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In peripheral heavy ion collisions the system has angular momentum. [T] It has been 
shown in hydrodynamical computations that the angular momentum leads to a large 
shear and vorticity [5]. Furthermore when the Quark-Gluon Plasma (QGP) is formed 
with low viscosity [3], interesting new phenomena may occur like rotation m , or 
turbulence, which shows up in form of a starting Kelvin-Helmholtz instability (KHI) 
016!. The deceleration of interpenetrating nuclei was observed and analyzed early 
in Ref. [7|. This leads to a rapid initial equilibration and to the development of a 
compact initial system. In peripheral collision this leads to considerable initial shear 
and vorticity, as well as to an almost complete conservation of the initial pre-collision 
angular momentum for the participants. 

Based on Refs. [210 we can extract some basic parameters of the rotation 
obtained with numerical fluid dynamical model PIGR. These parameters are extracted 
from model calculation of a Pb-|-Pb collision at = 2.76 TeV/nucl. and impact 

parameter b = 0.7bmax, with high resolution and thus small numerical viscosity. Thus, 
in this collision the KHI occurs and enhances rotation at intermediate times, because 
the turbulent rotation gains energy from the original shear flow. The turbulent 
rotations leads to a rotation profile where the rotation of the external regions lags 
behind the rotation of the internal zones. This is a typical growth of the KHI. 

The time dependence of some characteristic parameters of the fluid dynamical 
calculation were analysed in Ref. 0 . It was observed that R, the average transverse 
radius, Y, the longitudinal (rotation axis directed) length of the participant system, 
9, the polar angle of the rotation of the interior region of the system, are increasing 
with time. R and Y the speeds of expansion in transverse and axis directions are also 
increasing with time, while to the angular velocity of the internal region of the matter 
during the collision is decreasing. 

The initial angular momentum of the system is large, Ly = —1.05 x lO^S. As 
this is arising from the 2 ; directed beam velocity, initially at the vertical, x, edges the 
velocity difference is large, while horizontally the rotation starts up delayed, because 
this is not a solid body rotation. Here we considered the rotation measure versus the 
horizontal, z axis which starts up slower and reaches a maximum around 5 fm/c after 
the start of the fluid dynamical evolution, i.e. around 8 fm/c after the initial touch of 
the nuclear surfaces. 

Exact models, see e.g. Ref. 0, provide good insight into physical phenomena. 
We want to use the above mentioned fluid dynamical calculations to test a new family 
of exact rotation solutions of fireball hydrodynamics 0iini. This model offers a few 
possible variations, here we chose the version lA to test. We use the axis labeling 
of Refs. 0 [TU] . so that the the axis of the rotation is y while the transverse plane 
of the rotation is the [x, z] plane. Thus the values extracted from the results of the 
fluid dynamical model 0, should take this into account. The initial radius parameter, 
R, corresponds to the x axis in hydro, and we assume an a;, z symmetry in the exact 
model. The rotation axis is the y axis in hydro. The exact model assumes cylindrical 
symmetry, so it cannot describe the beam directed elongation of the system, but this 
is arising from the initial beam momentum, and we intend to describe the rotation of 
the interior part of the reaction plane and the rotation there. 

For simplicity we also assume that the Equation of State (EoS) is 
e = Kp and p = nT , 
with a constant k. 


( 1 ) 
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Now we calculate the equation of motion, (15) in Ref. [lOj. and its solution 
nm{dt + V ■ V)u = — Vp 
For the variables of this equation we have: 


( 2 ) 



n 



Vo , , 

noyi'is), 

1 -± I 


and in addition in Ref. m it is assumed that the temperature and the density have 
time independent distributions with respect to a scaling variable: 



If we asume cylindrical symmetry and use the coresponding cylindrical coordinates 
instead of {x,y,z), we can use the coordinates of length dimension, {rp,r^,ry), so 
that 


rp =p, ry, = rpp, Ty = y . 


These are the ’’out, side, long” directions. The characteristic values of these 
coordinates are then {R,S,Y). Then the scaling variables are introduced as 



where S is the roll-length on the outside circumference, starting from (/jq = 0 and 
So = 0 at to) S = Rip and ip = to and this displacement is orthogonal to the 
longitudinal and transverse displacements. The internal roll-length Vy, = iprp, the 
corresponding velocity is = tor p, and so r^. On the other hand from the 

scaling of rp, it follows that r^ = R^ Sp. 

In case of these scaling variables the distributions of density and temperature, n(s) 
and Tlt,s) should not depend on Sy, or r^p, just on the radius and the longitudinal 
coordinates. Therefore just as in Ref. [S] we introduce another scaling variable: 


S = Sp -I- s. 


Our reference frame is then spanned by the directions (rp,r,p,ry). In this case due to 
the cylindrical symmetry the derivatives, dsjdry, vanish. In this coordinate system 
the volume is P = ■nR?Y. 

Now following Ref. [5], for the right hand side of Eq. (|^: For the r.h.s. of this 
equation we have: 

— Vp = — VnT 



( 4 ) 
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TriV~^ 

where 7 = I/k and Q = ^ 


Using the p, ip, y coordinates, the rotation would show up as an independent 
orthogonal term. However, (as discussed in the Appendix) the closed system has no 
external torque, and the internal force from the gradient of the pressure is radial, 
which does not contribute to tangential acceleration. The change of the angular 
velocity arises from the angular momentum conservation in the closed system as a 
constraint, so we do not have to derive additional dynamical equations to describe the 
evolution of the rotation. 

Now for the left hand side of Eq. (i): the velocity field scales as 

V = VpBp-v^pe^+Vze^ = j^rpep-Lorpe^+yryBy . (5) 

We first calculate the time derivatives for the components. (See e.g. [11)1: 


dtVp 




— LJ 


r 




dtv^ 


— UJ 



dtVz 


Y 

Y 


Y2 

^72 


( 6 ) 


The other term of the comoving derivative includes 
^ d d d 

V -V = Vp- -h v^- -h VyT— 

dr„ ^ " dr„ 


and this term gives: 




R 


y 2 


{v ■ V)v = ^rpBp + uj-rpe^ + y^^VyBy 
By adding Eq. ([^ and (|^ we get: 


mn{dt+v ■'V)vp = mn (^R/R^ 
mn{dt+v • V)uy = mn (y/Y^ 1 


— UJ 


y ’ 


(7) 

( 8 ) 

(9) 


As a consequence the equality of the right hand side and left hand side of the 
Euler equation ([^ leads to the ordinary differential equations. Multiplying the two 
non-vanishing equations with R^ and Y^ respectively yields: 


RR - W/R^ = YY = 


Q 

(7ri?2y)7 ’ 


( 10 ) 


where W = i?g. From the angular momentum conservation uj = ojqRq/R^, and 
the rotational term, that appears in the equation, takes the form W/R?. 

Due to the EoS the pressure is proportional to the baryon density n, just as the 
r.h.s. of the Euler equation, therefore the equation of motion does not depend on n 
or ng. 
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Following Ref. [S], we want to calculate the total energy of the whole system, then 
we should integrate it for the whole volume, V. Thus, not only the scaling of 
V = {vp,Vip,Vz) but also the particle density distribution, n{s) will be considered. 


At the surface the rotational energy is Sside = \mS^ = and if we 

express uj via wq by the relation uj = ujqRq/R?, then Eside = WjR^, as before. The 
expansion energy at the surface is Eout = \rnR^, and for the longitudinal direction 
we have, Elouq = \rnY'^. 

We can calculate the radial and longitudinal expansion velocities and the 
corresponding kinetic energies, and also the kinetic energy of the rotation. In the 
evaluation of the internal and kinetic energies, the radial and longitudinal density 
profiles of the system are taken into account. 

Let us assume that the temperature profile is flat, and thus that the density 
profiles are Gaussian and separable. Further if we assume that the system size is 
finite, then the scaling variables, Sp and sy, may extend from 0 to 1. In this case at 
the external boundary we have to apply the necessary boundary conditions so that 
the solution of the Euler equation ([^ remains valid. With this approximation we 
calculated the different integrated energies (and shown in Appendices [Spl]) . 

Summing up the kinetic energies yields 

Ek^u = \mNB {a^R^ + a^RW + , ( 11 ) 

where in case of finite extent of the system 
a'^ = 4:y^CnlB{^SyM)Ic{^SpM) and 
P‘^ = 4:V2CnlA{lspM)lDilsyM) where 

Cn = l/[2V2lA{spM/2)lBisyM/2)], See Ref. [E] I (in terms of the integrals 
evaluated in Appendices [Spl] ). 

Here with {spM = {syM = 1, and /3^ are clearly time independent, because 
they depend on the scaling variable only, and we get the values = 0.4585 and 
^2 = 0.2911. 

Alternatively one can assume that the system size is infinite so that the scaling 
variables range from 0 to oo. In this case the Radius, and Length parameters, R 
and F, are considered as the width of the Gaussian scaling distribution. Thus the 
parameters will be 


Ia = Ib = 1, Ic = ‘^Id = a/tt 


( 12 ) 


and consequently c? = 2.0 and 0^ = 1.0. In the present case we follow this 
configuration. 

If we divide this result by the conserved baryon charge, Nb-, we will get 


dEkin 


-m 


a 




(13) 


Based on the EoS, e = Kp = nnT, we can calculate the compression energy based 
on the density profiles of n{s) and e{s) = Kn{s)T{s). 


t Ia{u) = 1 - exp(-n), Ib(u) = Ic(u) = 1 - (1 + li)exp(-M), Id{u) = ^#(07) - 

where $(ri) = erf(ii) = fg exp(—x^) dx. 
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Let us make the same simplifying assumptions on the density profiles as we did 
earlier. Now we will have the same density profile, normalized to Ng, for the volume 
integrated internal energy and the net baryon charge: 


Eint = K pdV = K nTdV 


= KNeToiVo/VrCn ^ [ Y ( v{s) dsp ^ 

V Jo Jo 


= kNbTo{Vo/V)'* = KiTiQ 
where Cn is the normalization constant. 


0 Jo 
1 

{TTR^Y)'r ’ 


(14) 


4. Reduction to a Single Differential Equation 


Following the method of Ref. m , we study the following combination of variables: 
E = {a^R^ + I3^Y^) 

= dt (a^RR + p^YY^ 

= a'^R^ + j3^Y^+ a^RR +P'^YY, (15) 


where we used the notation dt = % and 9^ = ^. We can replace the last two terms, 


a^RR, eYY, by using Eqs. (10), i.e. we use the Euler Eq. (^2|. Then we obtain: 


E = a^R^ + / 32 y 2 


2^ / 




Q 


(ttR^YP ’ 


(16) 


At the same time from the energy conservation, Efot = E/^tn + Eint, we get 

that 


Etot 

Ng m 


a^R^+P^Y^+a'^ 


W 

w 


2kQ 

JnWYfr 


(17) 


where we used the EoS and the parameter n now appears in the expression of the 
energy. If our EoS is such that 

K=^(a^+P^) = ^, (18) 

then E = 2Etot/{Ng m) = const., and in the same type of calculation as in Ref. |13) . 
we can introduce 


U^{t) = a^R^{t) + j3'^Y'^{t), 

which satisfies 


d'i {a^R^ + p'^Y'^) = d^U^it) = 2E . 

Thus, the solution of Eq. (20), we can be parameterized as: 


(19) 

( 20 ) 


U^{t) — “1“ — to) “1“ ^ 


( 21 ) 
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A = a^Rl + + a^W/Rl + 

B = 2a'^RoRo + 2/3‘^YoYo 


C = a^Rl + /J^yp • 

Due to the difficulties described in Appendix I 


( 22 ) 

we cannot use the method 

described in m- Instead let us take one of the Euler equations from Eq. ( |l0| ) , 

= YiSvfi • 

and express R? in terms of W^it) which is known based on the energy conservation; 

R^{t) = {U‘^{t) - ^'^Y'^)lc? , (24) 

and this will lead to the second order differential equation for Y (t): 


Y = 


which can be solved 


y [7rr(172(t)-^2y2)]' 


= fiYA) 


Then R{t) and R{t) are given by Eqs. (24) and (16) respectively. 


(25) 


t 

(fm/c) 

y 

(fm) 

y 

(c) 

R 

(fm) 

R 

(c) 

UJ 

(c/fm) 

0.0 

4.000 

0.400 

2.500 

0.250 

0.150 

1.0 

4.440 

0.469 

2.859 

0.704 

0.115 

2.0 

4.922 

0.490 

3.405 

0.834 

0.081 

3.0 

5.415 

0.495 

4.079 

0.877 

0.056 

4.0 

5.912 

0.497 

4.833 

0.894 

0.040 

5.0 

6.409 

0.497 

5.636 

0.902 

0.030 

6.0 

6.906 

0.498 

6.469 

0.906 

0.022 

7.0 

7.404 

0.498 

7.322 

0.909 

0.017 

8.0 

7.901 

0.498 

8.190 

0.911 

0.014 


Table 1. Time dependence of characteristic parameters of the exact fluid 
dynamical model [8]. R is the transverse radius, Y is the (rotation axis directed) 
length of the system, R, Y are the speed of expansion in transverse and axis 
directions, and lo is the angular velocity of the matter. 


The derivatives, i?(to) and y(to) in this exact model do not equal the ones 
obtained from the fluid dynamical model, because in the more realistic fluid 
dynamical model the density and velocity profiles do not agree with the exact 
model’s assumptions. Also initially in the realistic fluid dynamical model the angular 
momentum increases in the region due to the developing turbulence, while in the exact 
model the angular velocity is monotonously decreasing due to the scaling expansion. 

The Runge Kutta [14] method was used to solve this differential equation. We 
chose the constants, Q and W, as well as the initial conditions for R and Y. 

Based on the fluid dynamical model calculation results we chose the parameters: 
Tq = 250 MeV, m = 939.57MeV wp = 0.15 c/fm. For the internal region we take 
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Figure 1. (Color online) The time dependence of the kinetic energy of the 
expansion, the internal energy, Ei„t, the rotational energy, Erot, and the 

total energy, Etot per nucleon in the exact model with the initial conditions Rq 
= 3.5 fm, Yq = 5.0 fm, i?o = 0.25 c, Yg = 0.30 c, loq = 0.1 c/fm, k = 3/2, Tg 
= 400 MeV. For this configuration Etot = 816 MeV/nucl. The kinetic energy of 
the expansion is increasing, at the cost of the decreasing internal energy and the 
slower decreasing rotational energy. The rotational energy is decreasing to the 
half of the initial one in 3.3 fm/c. 


the initial radius parameters as Rq = 2.5fm and R = 0.25c, and we disregard the 
larger extension in the beam direction, because our model is cylindrically symmetric 
and because the beam directed large elongation is a consequence of the initial beam 
directed momentum excess. In this exact model the rotation axis, denoted by Y, 
corresponds to the out of plane, y direction in the fluid dynamical model (and not 
to the beam direction!). Due to the eccentricity at finite impact parameters, with an 
almond shape profile, the initial out of plane size is larger the in plane transverse size, 
so we chose initially Yg = 4.0fm and Y = 0.4c just as in Ref. [5]. (Table 

As the exact solution is able to describe the monotonic expansion, and so the 
steady decrease of the rotation, we start from a higher initial angular velocity than 
shown by the fluid dynamical model, PICK, as the angular velocity, measured versus 
the horizontal plane, starts from zero. 

Applying these initial parameters the exact model yields a dynamical development 
shown in TableAccording to expectations the radius, i?, and the axis directed size, 
Y, are increasing, the angular velocity, to decreases. The total energy is conserved, 
while the kinetic energy of expansion is increasing, and that of the rotation and internal 
energy are decreasing. See Fig. 
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The change of the expansion velocity, vr = R, is shown in Fig. [^left. The more 
rapid velocity change arises partly from the centrifugal acceleration of the rotation, 
but also from the fact that the initially smaller transverse size increases faster in the 
direction of equal sizes in both directions. See Fig. [fright. 

The study of the rotation in an infinite system is. on the other hand, problematic 
as we assume solid body rotation (i.e. the angular momentum applies to the whole 
infinite system). So the applicability of this infinite model to a heavy ion reaction is 
highly approximate, and the external tails should be disregarded. 

Other finite scaling expansion profiles can also be studied, based on the given 
examples, and these may fit more detailed fluid dynamical models better. 

5. The vorticity 

In the usual convention in heavy ion physics, the beam axis is the 2 ; axis, the impact 
parameter vector, b, points in the x direction, and the projectile is at positive x and 
moves in the positive 2 direction. Thus the rotation axis is the y axis, this is the axis 
of the cylindrical symmetry of the rotating exact model system we discussed above. 
The reaction plane x, z is spanned by the cylindrical coordinates r, (p in the discussion 
above. 

For rotation around the y-axis the vorticity is defined in terms of the velocity 
field as ojy = {dzVx — dxVz)l2. We use the conventions of the exact model here, so that 
we chose the rotation axis to be the y axis and the plane of the rotation is the x,z 
plane, which corresponds to the reaction plane. We assume that the rotating system 
is symmetric, so we introduce cylindrical coordinates around the rotation axis. 

For this configuration in cylindrical coordinates the vorticity is . 


u} = rotv 



t (fm/c) t (fm/c) 


Figure 2. (Color online) Left: The time dependence of the velocity of expansion 
in the transverse radial direction, vji and in the direction of the axis of the 
rotation, vy for the configuration shown in Fig. ^ The expansion velocity is 
increasing in both directions. While in the axis direction the velocity increases 
from 0.3 c to 0.6 c in 8 fm/c time, the radial expansion increases faster, in part 
due to the centrifugal force from the rotation. Right: The time dependence of the 
Radial, R, and axis directed, T, size of the expanding system. As the Y directed 
velocity is initially larger its change is relatively smaller. 
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n(Y,R) (c/fm) at t = 6fm/c n(Y,R) (c/fm) at t = 8fm/c 



Figure 3. (Color online) The energy weighted vorticity in the classical rotating 
exact model with Gaussian density profiles, with an EoS of k = 3/2, and with 
initial parameters: 3.5 fm mean radius, 5.0 fm mean length, 0.1 c/fm angular 
velocity, 0.25 c radial velocity, 0.3 c axis directed velocity (vy). The initial 
temperature of the matter is T = 400 MeV. The figures show the configuration 
at different times. At t = 2 fm/c, the mean radial (longitudinal) sizes and speeds 
are 4.27 fm (5.75 fm) and 0.494 c (0.437 c), and the angular velocity is 0.07 c/fm. 
The V = c-boundary is at Ymax = H-O fm and Rmax = 6.3 fm. At t = 4 fm/c, 
the mean radial (longitudinal) sizes and speeds are 5.67 fm (6.76 fm) and 0.62 
c (0.52 c), and u) = 0.04 c/fm. The v = c-boundary is at Ymax = 7.1 fm and 
Rmax = 11-1 fm. At t = 6 fm/c, the mean radial (longitudinal) sizes and speeds 
are 6.65 fm (7.91 fm) and 0.69 c (0.57 c), and uj = 0.0 c/fm. The v = c-boundary 
is at Ymax = 12.0 fm and Rmax = 7.98 fm. At £ = 8 fm/c, the mean radial 
(longitudinal) sizes and speeds are 8.04 fm (9.14 fm) and 0.73 c (0.60 c), and 
to = 0.02 c/fm. The v = c-boundary is at Ymax = 13.16 fm and Rmax = 9.32 fm. 


\p dip dy J ^ \dy dp ) \P dp p dp ) ^ 

= — 2uiey . ( 26 ) 

At the last step we use Eq. ([^, where Vy does not depend on p, Vy, does not depend 
on y, Vp does not depend on y, Vy does not depend on p and Vp does not depend on 
(/?, thus only one term contributes to the vorticity, which is directed in the direction 

of By. 

Thus the vorticity in this model is spatially homogeneous, and depends on the 
time only, uj = u}{t). However, from the point of view of observations, it is important 
what amount of energy or mass is representing a given fluid element with the given 
vorticity. In the solution we presented here we assumed a uniform temperature, which 
led to a gaussian density and energy profile, e(p) = KTn(p). 
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t (fm/c) 


Figure 4. (Color online) The time dependence of the kinetic energy of the 
expansion, the internal energy, Eint, the rotational energy, Erot, and the 

total energy, Etot par nucleon in the exact model with the initial conditions Rq 
= 2.5 fm, Vo = 4.0 fm, Rq = 0.20 c, Vq = 0.25 c, ljq = 0.1 c/fm, k, = 3/2, To 
= 300 MeV. For this configuration Etot = 576 MeV/nucl. The kinetic energy of 
the expansion is increasing, at the cost of the decreasing internal energy and the 
slower decreasing rotational energy. The rotational energy is decreasing to the 
half of the initial one in 2.9 fm/c. 


Following reference [5], we define an energy-density-weighted, average vorticity as 
^zx = w{rp, r^, Ty) UJ (27) 


so that this weighting does not change the average circulation of the layer, i.e., the 
sum of the average of the weights over all fluid elements is unity, {w{z,x)) = 1. This 
weighting does not change the average vorticity value of the set; just the cells will 
have larger weight with more energy content. 

Let us fist calculate the internal energy for a finite system: 

p+aY pbR 

E^rit = / e(rp,ry) 27vrpdrpdry , (28) 

J-aY Jo 


while the the energy at a given radius (at a given time) is e(rp,ry). Thus the weight 
density will be 


w{rp,r^,ry) 


T^^jrp.ry) 

EtotIV 


(29) 


where 


r°° = e [(1 + 1 /k) 7 ^ - 1 /k] and 
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(^irp,ry) = KUoyTo ly j T{s)e => ■'o 


= kToUq ( — 


Fn ^ 


g-Sa/ 2 g-Sp /2 


(30) 


while 7 ^ = 


1 - ^7 


(s^p) ~ (‘^’’p) ~ (F^y) ■ Here we assumed that T(s) = 1 

and let C = kuo^Tq to get e{rp,ry) = With s = Sp + Sy where 

Sp = rl/E? and Sy = 

We also get that 


= C 


/-\-<iY' ^2 phR 

e~2Y^dry / < 

-ay ./o 


e 2 ^ 2'Krpdrp . 


Using a change of variables to Sy and Sp, so that dvy = ^y^dsy and 2Trrpdrp = 
R^irdsp, the scaling integration boundaries will be SyM = SpM = and we find 


r+a-^ 


= CttR^Y 


„-sj2 


Sy Jo 


e ’"’‘/‘^dsp . 


We can express the integrals as follows 


|•b^/2 

2/a(&V2) = / e-’^o/^dsp = 2 

Jo Jo 


e '^du 


and 


where 


2y2/B(aV2)= r e-««/2^ = 2 T e-^»/2^=2V2 [ 

J—a^ \f^ Jo Jo 


lA{u) = l-e and /b(m) = v^$(v^) • 


■ 0^/2 


Thus 


= C^i?2F2/^ (62/2)2^2/3 (aV2) 

^F4V2 


= CttR 


a/v^ 1 . 


(32) 


(33) 


Similarly we can calculate the Kinetic energies, Ekim for the rotation and radial 
and longitudinal expansions as in [5] (and see the Appendix). Then, using Eq. 0 
or (17), the total energy of the system is 

^tot — ^int “1“ ^kin ? (34) 

which we can use in the calculation of the weighted vorticity. 

In the present study we assumed an infinite system with scaling gaussian density 
profile, Eq. (12), so that the integrals are evaluated up to infinity. 
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We performed a set of calculations to study the applicability of the model to heavy ion 
reactions. This presented non-relativistic model leads to super-luminous velocities at 
late times and at the external surface of the system. We used a parametrization where 
the peripheral energy density is cut off exponentially, and took initial conditions such 
that the vast majority of the system is in the non-relativistic applicable domain of the 
model. 

A first series of calculations is presented in Figs. 

With the parameters as defined in Figs. |^when we reach t = 8 fm/c the surface 
speed reaches the speed of light already when the Energy weighted vorticity drops to 
40% of the top central value. Thus, a substantial amount of matter is outside the 
range of physical applicability of the model. The evaluation of polarization would not 
be realistic with these sets of parameters. Therefore we modified our initial conditions 
such that the applicability of the non-relativistic model holds up to the final, freeze 
out time of about 8 fm/c. 

The time development of the change of the different forms of energy are presented 
in Fig. 1^ for the modified initial state. While the sizes R and Y, and the expansion 
velocities in these directions are shown in Fig. 

Now the total energy of the system is about 70% of the previous example, and 
the initial rotational energy is 60% of the previous one. 

We also performed another test series with this more compressed initial state 
configuration]^ While up to t = 6 fm/c the majority of the energy weighted vorticity 
is in the applicable domain (where the velocity does not exceed the velocity of light), 
at t = 8 fm/c roughly 95% of the energy content is still in the applicability domain 
of the non-relativistic exact model (see Fig. [^. We may estimate that about 50-70% 
of the initial energy of a peripheral collision will contribute to the expansion of an 
symmetric solution of our participant system. Thus the model is applicable at lower 
energies, FAIR and NICA, energies, while at the top energies of RHIC or LHC the 



Figure 5. (Color online) Left: The time dependence of the velocity of expansion 
in the transverse radial direction, and in the direction of the axis of the 
rotation, vy, for the configuration presented in Fig. The expansion velocity 
is increasing in both directions. While in the axis direction the velocity increases 
from 0.25 c to 0.5 c in 8 fm/c time, the radial expansion increases faster. Right: 
The time dependence of the Radial, R, and axis directed, Y, size of the expanding 
system. 
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0(Y,R) (c/fm) at t= 2fm/c fJ(Y,R) (c/fm) at t= 4fm/c 



Figure 6. (Color online) The energy weighted vorticity in the classical rotating 
exact model with Gaussian density profiles, with initial parameters as given in 
Fig. The figure shows the configuration a,t t = 2 fm/c, when the mean radial 
(longitudinal) sizes and speeds are 3.16 fm (4.64 fm) and 0.43 c (0.38 c), and the 
angular velocity is 0.06 c/fm. The boundary is at the position where the velocity 
of matter reaches the speed of light, c. This happens at Ymax > 9.0 fm and 
Rmax = 5.68 fm. At t = 4 fm/c, the mean radial (longitudinal) sizes and speeds 
are 4.11 fm (5.51 fm) and 0.54 c (0.45 c), and oj = 0.04 c/fm. The v = c-boundary 
is at Ymax = 10.58 fm and Rmax = 6.24 fm. At t = 6 fm/c, the mean radial 
(longitudinal) sizes and speeds are 5.25 fm (6.45 fm) and 0.59 c (0.48 c), and 
u) = 0.02 c/fm. The v = c-boundary is at Ymax = 11.34 fm and Rmax = 7.56 
fm. At t = 8 fm/c, the mean radial (longitudinal) sizes and speeds are 6.44 fm 
(7.49 fm) and 0.62 c (0.51 c), and uj = 0.02 c/fm. The v = c-boundary is at 
Ymax = 12.58 fm and Rmax = 8.75 fm. 


reliability of this model is qualitative, and may provide estimates with 15 - 20 % 
accuracy. 

7. Conclusions 

The effect of QGP formation on the directed flow and the arising 3rd flow component 
or antiflow was first observed in fluid dynamical calculation at energies above 10 GeV 
per nucleon in Ref. m- The nuclear EoS has to satisfy strong constraints from 
the observed Neutron and Hybrid Star masses [H] Spin-orbit interaction and the 
momentum dependence of the nuclear interaction m, influence the nuclear EoS and 
developing rotation and polarization of the participant matter. The nuclear EoS has 
a strong effect on the collective motion. Transverse flow and collectivity was observed 
early both in fluid dynamical, nuclear cascade and molecular dynamics models m- 
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In conclusion, the exact model can be well realized with parameters extracted 
from detailed, high resolution, 3+ID relativistic fluid dynamical model calculations 
with the PICK code. It provides an estimate of the rate of decrease of angular speed 
and rotational energy due to the expansion in an explosively expanding system. This 
indicates that the effects of rotation can be observable in case of rapid freeze out and 
hadronization, although the Kelvin Helmholtz Instability is not present in this model 
and this reduces the rotation at later times. 

This indicates that the presence of the KHI is essential for an observable effect 
of the rotation, and thus the observation of the rotation is strongly connected to the 
evolving turbulent instability in low viscosity Quark-gluon plasma. 
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Let us evaluate the baryon density, n{s), and for simplicity let us assume that in 
case lA of Ref. m the temperature is constant, T(s) = 1, then it follows that, 
i^{s) = {Nb/V) exp(— s/ 2), where Nb = no Vq- Due to the exponential density profile, 
if s is a sum of the coordinates in two orthogonal directions, as s = Sp = Sz, then n(s) 
separates into two multiplicative terms: np(sp) and nz(sz). For further simplifying 
the formalism, we can introduce a coordinate change s = 2u for integrals of type 
lo^"^ /(u/2) du. Then ds = 2du and 2 fin/2) du = /(s) ds. This change will 

thus modify the upper limits of integration, and the normalization by a factor of two. 
These adjustments are included in the final expressions in Appendices [8pT] 

The baryon density distribution is then 


( \ AT IY'^ 

n{rp,ry) = Nb-^b e y' 


where Cn is a normalization constant, which will be determined later. The 
normalization can be performed up to a finite size, i? and Y, or up to infinity. 

Cr,. 


27r 


irp,ry) = Nby- 

pB./oo 

( 

f'Y/oo 


/O 


-^M^\pdTp 


' —Y/oo 


e dry , 


(35) 


here the first integral up to infinity gives r(l)i?^/2, while the second one r(0.5)y = 
a/ttY . In X, y, z coordinates this is: 


) = Nb X const. 




Or in cylindrical coordinates 


Nb^x 7rr(l)R2 X = Nb^tt^^^R^Y , 


(36) 


(37) 


which is the same. The integrals were evaluated up to limits in infinity. If we perform 
the definite integrals up to a finite limit, we get similar scaling behaviour. Let us now 
change the variables to scaling variables introduced in Ref. m, but in cylindrical 
coordinates. 

n(sp, Sy) = 

Now using the relations 2Trrpdrp = irR^dsp and dry = ^^^ dsy we get 


(sp, Sy) — Nb y 


dsu 


X ttR^ / e dspXY e ^ = 




7rR2y 

' V 


/■U.t) rU.t, , 

X 2/ e-^duxV2 e-" — 

. n . n \/u 
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= NbCJ^ 2 /^(0.5) V2Jb( 0.5) 
This should be equal to Nb thus the normalization constant is 


Cn = l 


2v^/a(0.5)/b(0.5) 


(38) 


Here /^(0.5) and /b( 0.5) are constants, which do not change during the scaling 
evolution, when the density profile remains the same. At infinity /^(oo) = r(l) 
while /b(cx)) = r(0.5) = but at different integration limits the ratio of the two 
integrals will be different m- 

Ia(,u) = 1 ~ exp(—n) 

/b(u) = V^$(V^) , (39) 


where 

2 /■“ 

$(ri) = erf(ri) = / exp(—a;^) da; . (40) 

Jo 


9. Appendix - The Moment of Inertia 


Consider a body with scaling expansion, and with solid body rotation (i.e. the angular 
velocity is uniform for the whole body, uj = uj{t) but it does not depend on the spatial 
coordinates. Let us denote the moment of inertia with 0, 

0 = y mn{r)r^dJr. (41) 

Then the angular momentum and the rotational energy are 

L = Qui and E = 

Now we assume that our system has no external torque, and all internal forces are 
radial, so the angular momentum must be conserved, during the scaling expansion 
driven by the pressure gradient which is radial in a cylindrically symmetric system. 
Thus, the angular velocity is not directly influenced by the dynamics, just via the 
angular momentum conservation. From L = 0 , it follows that 


0 Cli = —0 w or u) = ■ 

Thus the change of the angular velocity is a direct consequence of the change of the 
moment of inertia 0, while 0 is proportional with the square of the radius of the 
system in a scaling expansion where the density profile remains the same during the 
expansion. Consequently 


0 

= = ““TP 


UJ — UJQ- 


Ro 

Id 


We still have to evaluate the moment of inertia accurately to provide precisely the 
energy of rotation. Thus, using the scaling variables 


0 = 

V 


3 X / e 


dsv 


= mNBR^ Cn 4 Jc( 0 . 5 ) 72/5(0.5) , 


( 43 ) 
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where Ic{u) = 1 — (1 + u) exp(— m) . 

As before these integrals do not change during the scaling expansion, on the other 
hand the volume and the moment of inertia have different coefficients in the energy 
expression. As a consequence the kinetic energy of the rotation is 

Emu = = \mNB C„4y2/c(0.5)/B(0.5)i?2a;2 

= ia^TOAfBi?V. (44) 

Here we have introduced the constant 

a^ = ^y/2CM\)Ici\), (45) 

that can be used in the main course of the work. 


10. Appendix - Kinetic energy of radial expansion 


The radial velocity is given by Vp = {R/R)rp and consequently = R^Sp. Thus the 
kinetic energy of radial expansion is 


tuNb ttR^Y -2 

-^kin — 2 ^ 

tuNb 


[ e ^p^^Spdsp [ e 

/o Jo 


C„4y2/c(0.5)/B(0.5) 


R^ 


(46) 


11. Appendix - Kinetic energy of longitudinal expansion 


The longitudinal velocity is given by Vy = (Y/Y)ry and consequently Vy = 
Thus the kinetic energy of longitudinal expansion is 

Ekin 

tuNb 


Cn 4^2/^(0.5) Jc (0.5) ^ 


where Id{u) = — ^/ue Here we can introduce the constant 

which will be used in the calculation. 


Y^Sy. 


(47) 

(48) 


12. Appendix - Not realizable analytic solution. 


One may find a solution for the dynamical evolution of R{t) and Y(t) based on Eq. 16 


with simplifying the problem to a singe, first order differential equation in a similar 
way as it is done in Ref. m- Let us introduce a parametric function, ()>(t), 


aR{t) = U(t) sin (j){t) 
(3Y{t) = U{t) cos 4>(t) , 


(49) 
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satisfying Eq. (19). Now inserting Eqs. (49) into Eq. (16), and noticing that 
= IP + we get the following first order differential equation for (j): 


1 

U^t)[ 


T- U^it) - 


a'^W 


U^(t) sin^^ (ttU^ (t) sin'^ (/) cos _ 


.(50) 


The initial value of the variable (j) is chosen such that Eq. (49) is satisfied for 
U{to), 

The problem with this solution is that Eq. (50) describes the square of (p and 
in a realistic situation it is not trivial to find the sign of the r.h.s. of the dynamical 
equation for (j). This sign alternates. 












